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Abstract 

A number of demographic factors, many of which related to human-driven 
encroachments, are predicted to decrease the effective population size (N e ) rel- 
ative to the census population size (AT), but these have been little investigated. 
Yet, it is necessary to know which factors most strongly impact N e , and how to 
mitigate these effects through sound management actions. In this study, we use 
parentage analysis of a stream-living brown trout (Salmo trutta) population to 
quantify the effect of between-individual variance in reproductive success on 
the effective number of breeders (N\,) relative to the census number of breeders 
(Ni). Comprehensive estimates of the N^/N ratio were reduced to 0.16-0.28, 
almost entirely due to larger than binomial variance in family size. We used 
computer simulations, based on empirical estimates of age-specific survival and 
fecundity rates, to assess the effect of repeat spawning (iteroparity) on N e and 
found that the variance in lifetime reproductive success was substantially higher 
for repeat spawners. Random family-specific survival, on the other hand, acts 
to buffer these effects. We discuss the implications of these findings for the 
management of small populations, where maintaining high and stable levels of 
N e is crucial to extenuate inbreeding and protect genetic variability. 



Introduction 

The effective population size (N e ) is arguably one of the 
most important parameters in evolutionary and conserva- 
tion biology as it determines the rate of genetic drift, the 
rate of loss of diversity of neutral alleles and the rate of 
increase of inbreeding experienced by an isolated popula- 
tion (Charlesworth and Willis 2009). This is because it 
establishes the relative evolutionary importance of sto- 
chastic (random genetic drift) and directional (migration 
and selection) factors and is thus closely being linked to 
genetic diversity (Crow and Kimura 1970). Genetic varia- 
tion is the raw material for evolutionary change, allowing 
populations to evolve as a response to various environ- 
mental challenges (Frankham 1996), be they natural or 



anthropogenic, and its maintenance is therefore a funda- 
mental concern in conservation biology. To be able to 
properly manage wild resources, genetic and evolutionary 
considerations have to be integrated into management 
practices (Young 2004; Waples and Naish 2009). 

Because the magnitude of drift is inversely proportional 
to the effective size (e.g., Crow and Kimura 1970), small 
natural populations or populations experiencing fragmen- 
tation and loss of connectivity, often the focus of conser- 
vation programmes, are vulnerable to the random loss of 
genetic variation. Knowledge of the effective size is there- 
fore particularly important in such cases (Lande and Bar- 
rowclough 1987; Nunney and Elam 1994). Anthropogenic 
activities have recently resulted in many wild populations 
being seriously reduced in size, leading to significant 
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reductions in N e (Frankham et al. 2010). Supportive 
breeding programmes may also be disadvantageous from 
a conservation genetics perspective if the individuals 
released into the wild originate from captive populations 
with small N e values (Ryman and Laikre 1991). Also, the 
efficiency of natural selection is directly proportional to 
N e because genetic drift is weaker when N e is large, and 
this will determine the fate of natural populations 
(Charlesworth 2009). 

Despite its fundamental importance in conservation 
genetics, it has been notoriously difficult to obtain reliable 
N e estimates in the wild. Accurate ecological N e estima- 
tion requires quantification of a number of demographic 
parameters that are difficult to obtain from natural popu- 
lations in most species. Examples of these include 
accounting for spatial substructuring and overlapping 
generations. In addition, a variety of factors related to its 
interpretation complicate its application in wild popula- 
tions. The ratio between the effective population size to 
the number of sexually mature adults in a population 
(census size), the N e /N ratio, is of great potential interest. 
Different groups of organisms, however, seem to be char- 
acterized by different NJN ratios, leading to confusion 
regarding the magnitude and the meaning of this parame- 
ter (Nunney 1993; Frankham 1995). Also, if there is a 
relationship between N and N e , variable population size 
can have a profound influence on long-term N e because 
effective size for a series of generations is approximately 
the harmonic mean of the single generation N e (Wright 
1938). This effect would be particularly relevant in popu- 
lations of conservational concern, where an apparently 
healthy census population may mask small N e (Nunney 
and Elam 1994; Frankham 1995). 

The most important factors responsible for reducing N e 
below the number of sexually mature adults in the popu- 
lation is a larger than binomial variance in individual 
reproductive success, unequal sex ratio and temporal fluc- 
tuations in population size (Wright 1938). Polygamous 
mating systems would result in high variance in repro- 
ductive success that in turn leads to low effective popula- 
tion size (Nunney 1993; Ardren and Kapuscinski 2003). 
In real populations, the variance of progeny number per 
family generally exceeds the mean progeny number in 
both sexes (Crow and Kimura 1970), so that the effective 
number of individuals of either sex is less than the actual 
number present in any generation. However, as genera- 
tions overlap, there is no easy way of quantifying this 
effect on the per-generation N e (Waples 2010). Estimates 
of reproductive success are commonly obtained using sea- 
sonal measures, but these primarily estimate the effective 
number of breeders per year (Nb) rather than the effective 
population size per generation (N e ). The Nb values are 
however only useful in estimating N e if it can be assumed 



that female and male fecundity is age-independent (Nun- 
ney 1993). Because N e is based on lifetime reproductive 
success, such estimates are much more difficult to obtain 
for wild populations of iteroparous organisms. Moreover, 
to estimate generational N e from demographic data, par- 
ents and their offspring should be counted at the same 
life stage. Family-correlated survival is expected to reduce 
Ne further compared to random survival (Crow and Mor- 
ton 1955; Hedgecock 1994). Little is currently known 
about the relationships between life history plasticity, 
demographic change and N e (Lande and Barrowclough 
1987; Nunney 1993), partly due to scarcity of the ecologi- 
cal data required for these analyses. 

Salmonid fishes exhibit a wide diversity of breeding 
systems (Fleming 1998; Fleming and Reynolds 2004), and 
even different populations of the same species living 
under different environmental conditions might have 
evolved different life history tactics and can thereby exhi- 
bit different breeding behaviours (Gross 1991). Both sexes 
in salmonids show high variance in reproductive success 
(Fleming 1998; Garant et al. 2001; Serbezov et al. 2010). 
However, reduced variance in individual reproductive 
success at low breeder abundance (genetic compensation) 
can lead to an increase in N e and might be a realistic 
aspect of salmonid breeding systems (Ardren and Kapu- 
scinski 2003; Frasier et al. 2007). Understanding the rela- 
tionship between N e and the mating system is therefore 
important in allowing comparisons between populations 
in conservation contexts (Nunney 1993). Proper manage- 
ment of such populations should account for those fac- 
tors that may potentially strongly impact (i.e., reduce) 
the N e and thus genetic variance. Indeed, in addition to 
climate change, salmonid fishes are exposed to a large 
number of human encroachments such as harvesting and 
fragmentation leading to reduced connectivity. Many of 
these encroachments influence the breeding system, levels 
of gene flow, population structure and overall fitness, 
thus impacting N e and the overall level of genetic varia- 
tion. To manage such populations, it is necessary to 
know which factors most strongly impact N e and how to 
mitigate these impacts through management actions. In 
this study, we use extensive demographic and genetic 
data on a brown trout (Salmo trutta) population to 
quantify the effects of variation in reproductive biology 
on Nb and N e . We use empirical estimates of age-specific 
survival (Carlson et al. 2008) and age- and sex-specific 
fertility (Serbezov et al. 2010) in an individual-based sim- 
ulation model to estimate the variance in lifetime repro- 
ductive success and obtain a generational N e estimate. 
Having this information, it is possible to evaluate what 
kind of population structure best maintains high and sta- 
ble levels of N e and thus help make informed manage- 
ment decisions. 
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Materials and methods 

Background 

The study site is a small forest stream, Bellbekken, in 
South-East Norway (N: 61°15', E: 11°51') that was sam- 
pled extensively during the period 2002-2007 (see Olsen 
and Vollestad 2003; Carlson et al. 2008; Serbezov et al. 
2010 for detailed information). Compared to nearby riv- 
ers, Bellbekken has limiting growth conditions for fish 
and is seldom, if ever, visited by anglers. The sampled 
section starts at a small waterfall, which should prevent 
upstream migration under most conditions, and there is a 
weak but significant genetic differentiation between trout 
in the stream and the larger river below the waterfall 
(Taugb0l 2008). 

The analyses described herein utilize empirical data 
described by Serbezov et al. (2010). Briefly, these were 
based on a total of 451 (195 males and 256 females) 
mature brown trout that were sampled (nondestructively) 
from Bellbekken during the 2002, 2003 and 2004 spawn- 
ing seasons, and 1856 juvenile trout sampled the follow- 
ing years, aged 0+ (i.e., during their first year of life) to 
2+ years and representing birth years or cohorts 2003 
(834 individuals), 2004 (707) and 2005 (315). All fish 
were genotyped at 15 microsatellite loci. The genotypes 
were used to individually assign (with 90% confidence or 
higher) offspring individuals to potential parents, using 
the MasterBayes software (Hadfield et al. 2006). Alto- 
gether 935 of the 1856 sampled offspring were assigned a 
father and 723 were assigned a mother (Serbezov et al. 
2010; table 3). Of the 451 potential parents, 78 males (of 
195) and 116 females (of 256) were assigned one or more 
offspring to them, whereas 117 males and 140 females 
were not assigned with any of the sampled offspring. We 
combined parentage assignment and sib-ship reconstruc- 
tion to improve our pedigree inferences as follows. The 
software COLONY v2.0 (Wang 2004; Jones and Wang 
2010) was used to partition the offspring cohort into full 
and half-sib families and to infer their parental genotypes. 
For the 2002 breeding season, a few large offspring half- 
sib families had no matching parents among the sampled 
adults (Serbezov et al. 2010). These were assigned to gen- 
erated parental genotypes and included in the calculations 
of the family-size variance as paternal half-sib families as 
they represent a large proportion of sampled offspring 
(Serbezov et al. 2010). 

Bellbekken has poor growth conditions for brown trout 
and age determination based on scale reading is therefore 
uncertain. Genotype matching was used to quantify the 
error in age determination for adult fishes. Fish some- 
times loose their Passive Integrated Transponder (PIT)- 
tags (Prentice & Flagg 1990), some of these individuals 
were tagged and genotyped again, and identical geno- 



types, combined with compatible length measurements, 
helped us confirm that an individual had been captured 
before. Specifically, we found that age was underestimated 
by 1 year in 35% and by 2 years in 18% of individuals 
that were 3 years or older. This uncertainty was incorpo- 
rated when calculating generation times and modelling 
maturation probabilities, as described below. Ageing of 
juveniles was based on length and is thus not affected by 
errors in scale readings. 

Statistical analyses 

We calculated the mean (fc) and variances (V k ) of off- 
spring numbers per male and per female parent for each 
spawning season using data presented in Fig. 1. Because 
offspring were observed and counted when juvenile (age 
0+ to 2+), and not at the same age as their parents, 
effects of mortality until reproductive age are not fully 
accounted for in these estimates of family size and their 
variances. To investigate the effect of family-specific mor- 
tality, we used the approach of Crow and Morton (1955) 
and adjusted k at different stages to their expected value 
of k = 2 (for a stable population) and recalculated V k . 
We calculated the mean (k) and variance (V k ) of the 
number of progeny for each sex for 1 -year-old offspring 
and for 3-year-old offspring separately. The age-specific 
ratios of variance to mean family size R t (= V kl /ki) and 
^3 (= ^W^s) values were then scaled to their expected 
values at k = 2 by assuming random survival from the 
last measured stage (Crow and Morton 1955; Waples 
2002), and these are denoted as R*i and R* 3 . This proce- 
dure yields an indication of how much family-size vari- 
ance through a particular life stage has contributed to 
the reduction of the effective population size (using 
N e /N = kl{\ + £* age ) (Waples 2002). 

Calculation of N e per generation requires an estimate 
of the variance in lifetime reproductive success (Hill 
1972, 1979). As we only have estimates of reproductive 
success for single seasons, we used an individual-based 
simulation model for a population of 98 males and 108 
females (the estimated average number of spawners in a 
spawning season for the 2002-2007 seasons, see below). 
We estimated maturation probability, m(a), at a given 
age, a, from the observed proportions, o(a), at that age 
(Barot et al. 2004): 



Uncertainties in observed age readings were included in 
the model by randomly sampling without replacement 
fish from the distribution of observed individuals and 
adding 1 or 2 years to the recorded age according to the 
observed error probabilities. This bootstrap technique was 
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Figure 1 Observed distribution of reproductive success (number of progeny assigned) for male and female brown trout in Bellbekken during 
three consecutive spawning seasons. 



performed 1000 times. Because we found only marginal 
differences between estimates for males and females, we 
pooled data for both sexes (Fig. SI). The individual-based 
simulations were done by following individuals through 
their reproductive life, from 3 to 8 years of age, as all 
assigned parents were within these ages. Yearly reproduc- 
tive success for mature fish was modelled by randomly 
sampling with replacement from the observed sex- and 
age-specific distributions of variation in reproductive suc- 
cess (Fig. 1). The number of offspring produced each year 
was then summed to give an individual estimate of life- 
time reproductive success. Survival was modelled using 
the age-specific survival rate previously estimated based 
on capture-recapture estimates from our system (Carlson 
et al. 2008). This bootstrap technique was also performed 
1000 times. 

Effective population sizes 

We estimated effective population sizes per season (N h : 
effective number of breeders) and per generation (N e ) 
from the seasons and lifetime distributions of family sizes, 



respectively. The effective numbers of breeders were cal- 
culated for each season and sex from the estimated mean 
number of offspring (k) and the corresponding variance 
(Vk) (Table 1). For females, this was done using: 

Nbfp) = W(i) ~ 1) , . 

(1 + Vkf/fe) u 

where N fi is the estimated number of sexually mature 
females in a spawning season i, and k ( and V kf are the 
mean and the variance of the number of progeny pro- 
duced (Crow and Denniston 1988). The equation for 
male spawners is analogous. The number of sexually 
mature males and females in any given year was estimated 
using the results from the three-pass electrofishing results 
and the standard Zippin depletion method (Zippin 1958; 
Bohlin et al. 1989). To obtain confidence intervals for the 
estimated Nb(i)'s, we resampled the observed number of 
spawners, with replacement, from their reproductive suc- 
cess. This bootstrap procedure was repeated 1000 times, 
and the 2.5% and 97% percentiles were taken as the 
lower and upper 95% confidence limits for Ntyi), respec- 
tively. 
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Table 1. Estimates of W b for 2002-2004 in Bellbekken based on sex ratio and variance in the reproductive success (based on parentage assign- 
ment analysis). Shown are the estimated numbers of males and females, Wb is the number of breeders corrected for unequal sex ratio (after eqn 
(1)). ^m, V m and kf, Vj are the mean and the variance in reproductive success for males and females, respectively, and Wb m and Nbf (o are the 
respective effective numbers of males and females each season (after eqns 2 and 3). 
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An estimate of effective number of breeders in a year 
was obtained by combining female and male effective 
numbers for that year: 



b(i) 



4(N bf(i) x N bm(i) 
(N bf(i) +N bm(i) ) 



(3) 



We also used the actual number of spawners (N m and 
Nf) in place of effective numbers in this eqn (3) to quan- 
tify the relative effects on N e of unequal sex ratio alone. 

The per-generation effective size, N e , was estimated 
from lifetime variances in family sizes obtained using the 
simulations described above and based on Hill's (1972, 
1979) approach, which is appropriate for species with 
overlapping generations. We use Hill (1979), eqn 10: 



N e = &NL/(a 2 m + a\ + 4) 



(4) 



where N is the average estimated number of spawners 
(males plus females) in a season, L is the average genera- 
tion length for males and females, and a 2 m and a 2 are the 
lifetime variances in reproductive successes for males and 
females, respectively. As presented by Hill (1972, 1979), 
calculations of these variances include separate variance 
and covariance terms for parent-offspring sex. As second- 
ary sex characteristics are not developed until the fish 
approach maturity, offspring were not sexed in our study 
and we substituted a 2 m and a 2 with Vj^, and V k[ , respec- 
tively. 



Cost of reproduction 

We investigated the effect of previous reproductive effort 
on later reproductive success (i.e., cost of reproduction) 
from estimates of individual reproductive success over the 
three consecutive spawning seasons. The individual repro- 
ductive success in a given season may depend on sex, 
body length, age, identity of the spawner and year 
(because of e.g. environmental differences among years). 
Body length and age are highly correlated, and we used 
body length as covariate in the analysis because we have 
higher confidence in the precision of body length mea- 
surements than age (cf. ageing errors above). We use a 
generalized linear mixed effects model implemented in 
the R package MCMCglmm (Hadfield 2010) with the fol- 
lowing structure to model effects on reproductive success: 

y = Pi Season + f> 2 Sex + ^Length + z^Spawner 
+ z 2 Year+e. (5) 

where Season (the number of the consecutive season an 
individual reproduces; from 1 to 3), Sex and fork Length 
are fixed effects, and Spawner identity and Year (2002, 
2003, 2004) are random effects, and £ is random residual 
error. The response variable (y), specified in the model 
with Poisson distribution, is the estimated number of 
progeny produced per individual and year. Significant 
negative values would be indicative of incurred costs 
of previous reproduction on current reproductive output. 



Sensitivity analysis 

To assess how the N e estimates depend on the values of 
the model parameters, we re-analysed the model with 
varying parameter values. First, we ran the simulations 
using the age-maturation curve (Fig. SI) unadjusted for 
ageing errors. Then, we ran the simulations with a range 
of values for the male and female variance in reproductive 
success (0.5-10 times the observed value) and age-specific 
survival probabilities (from 40% less to twice the 
observed value). 



Results 

Spawning brown trout found in Bellbekken were between 
2 and 8 years of age. There was a tendency for females to 
mature slightly younger than males (mean age 5.8 vs 
6.1 years; f 52 i = 3.9, p = 0.0006). The generation length 
of the two sexes, defined as their average ages at the time 
the offspring are born (the spring after a spawning event), 
was calculated from estimates obtained from the parent- 
age data. The estimated generation time for males was 
6.5 years (SD = 1.1 years) and for females 5.9 years 
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(SD = 1.0), for an average of L = 6.2 years (SD = 1.1), in 
good agreements with estimates from other Scandinavian 
stream-resident brown trout populations (I = 6.2 - 
6.6 years, Palm et al. 2003). 

Effective spawning population size 

The estimated census number of breeders (N) in Bellbek- 
ken from 2002 to 2007 was relatively stable, ranging from 
66 to 98 males and from 84 to 132 females. Based on 
these estimates, the effect on unequal sex ratios on N b for 
2002, 2003 and 2004 seasons was very small (eqn 3), and 
the iVb(sex ratk>)/N ratios were negligibly different from 
unity in all three spawning seasons (Table 1). 

In contrast, including the variances in family sizes (V k ) 
from the parentage assignment analysis yielded (eqs 2 and 
3) markedly reduced N b estimates, with N b /N ratios in the 
range from 0.16 to 0.28 (Table 1, right). The impact of var- 
iance in reproductive success can be considered separately 
for each sex and spawning season by comparing effective 
numbers (N\,{ and N bm ) with census numbers (Nf and N m ). 
From such comparisons, it is clear that deviations from 
binomial family-size distribution has a dramatic effect on 
reducing effective numbers for either sex (cf. Table 1), and 
the low N>/N ratios were not caused by a single sex only. 
This was observed in all three seasons and thus appears to 
be a general feature in this brown trout population. 

Effects of among-family variation in survival on esti- 
mated lifetime means and variances in family sizes were 
assessed by the ratio of variance to mean progeny num- 
bers (-R ; : Table 2), calculated separately for age 1 and 3 
progeny. R l values were always larger than 1, implying 
that family-size variances were larger than the binomial. 
This was more pronounced for males which, in addition 
to family-specific survival in the early juvenile stages, also 
reflect the reproductive skew in this population (see Serb- 
ezov et al. 2010). The adjusted R* t values (indicating 
what the index of variability would be if survival was ran- 
dom from that point on to an adult population with 
k = 2) were in general smaller than the unadjusted Ri 



values, but still considerably larger than unity. The com- 
bined effects of all demographic processes, including 
among-family variation in juvenile survival, led to a fur- 
ther reduction (in most cases) of Nb relative to those pre- 
sented in Table 1 and yielding final N b /N ratios between 
0.14 and 0.41 (Table 2). The reduction in N b /N between 
the 1+ and 3+ stages was smaller relative to the reduction 
up to the 1+ stage, ranging from 32.41% further reduc- 
tion (for females in 2002) to a 9.36% increase (females 
2003). The mean family size at enumeration at the 3+ 
stage (between 0.4 and 1.4; Table 2) is less than what 
would be expected in a population of constant size, 
k = 2. The Crow and Morton (1955) approach for scaling 
k and Vk> however, also applies if k < 2, the result being 
a scaled variance effective size that is equivalent to the 
inbreeding effective size that is calculated with the 
unsealed k and V k . All in all, the mortality between 1+ 
and 3+, although variable between sexes and seasons, 
does seem to be rather random among families, acting to 
reduce the negative effect of the strong skew in individual 
reproductive success on effective size and N b /N ratios. 

Repeat spawning 

A third of the spawners in Bellbekken were observed to 
be sexually mature during more than one season; a few 
individuals of each sex were even observed to be mature 
during five consecutive seasons (Fig. 2A). There was no 
difference between the sexes in the proportion of iterop- 
arous (repeat spawning) individuals (% 2 = 0.68, df = 1, 
p = 0.41) or in the levels of realized iteroparity (being 
assigned offspring in more than one season) (x 2 = 0.94, 
df = 1, p = 0.33) (Fig. 2B). On average, iteroparous males 
increased their reproductive success by 52% and females 
by 38% relative to semelparous fish. The only parameter 
from the mixed linear effect model (eqn 5) that had sig- 
nificant effect on reproductive success was body length 
(mean + SD: 0.09 ± 0.02). There seem to be quite large 
among-year variation in the number of offspring 
(mean + SD: 21.20 + 58.49) that masks any possible 



Table 2. Demographic estimates used in analysing family-correlated survival. Shown are the estimated census spawner numbers (N), mean (k) 
and variance (14) in family sizes for each sex and the corresponding age-specific index of variability (R-, = Vy^) and the scaled to their expected 
value at k = 2 age-specific index of variability (R*). These values are given for 2002 and 2003 and for the 1+ and 3+ age classes. The rightmost 
column shows the percentage change because of family-size variation that has occurred in the N b /N ratio between 1+ and 3+ stages. 
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Figure 2 (A) Proportions of brown trout individuals that were 
observed to have been sexually mature during different consecutive 
breeding seasons (1-5 seasons). The four categories (2002-2005) 
denote the first year a spawner was observed to be sexually mature. 
(B) Proportions of assigned spawners of both sexes in 1-3 consecutive 
seasons. 



effects of previous reproduction on the current offspring 
production, and no cost of reproduction was detected. 

Per-generation effective size 

Using the lifetime variances in estimated family sizes 
(<i 2 m = 82.6 and a\ = 22.7), L = 6.2 (the average genera- 



tion length) and N = 196 (average estimated number of 
spawning individuals for the 2002-2007 period) in Hill's 
eqn (4), N eV was estimated to be 104.3 ± 44.9 (SD), and 
a N eV /N ratio of 0.53 ± 0.24 (SD). 

To evaluate the consequences of repeat spawning on 
effective size in brown trout, we introduced obligate mor- 
tality (semelparity) after a maturation event in the model. 
This had the effect of reducing the number of offspring 
per breeder substantially for both sexes (3.8 + 8.8 vs 
1.8 ± 5.2 for males, f-test: p < 0.0001, 2.7 + 4.7 vs 
1.8 ± 3.7 for females, f-test: p < 0.0001), but also reduced 
its variance (F-test to compare two mean-standardized 
variances: ^1,13460 = 

2.8, p < 0.0001 for males; 
■Fi, 13480 = 1.6, p < 0.0001 for females). This reduction in 
variance in family size resulted in a N eV estimate for the 
hypothetical semelparous situation of 272.0 (SE = 112.7), 
more than twice the estimate with repeat spawning 
(Fig. 3). The increased family variance however led to a 
increased variance (uncertainty) in the N eV estimate (F- 
test to compare two mean-standardized variances: 
F li9 99 = 5.7, p < 0.0001). 

Sensitivity analysis 

The N eV estimates were relatively sensitive to changes in 
the variance in reproductive success, especially at lower 
values compared to the obtained values (Fig. 4A), but 
very little with respect to the value of the age-specific 
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Figure 3 Distribution of simulation-based estimates of A/ eV for the iteroparous model (A) and the semelparous model (B). The lifetime variance in 
reproductive success is simulated based on sex- and age-specific fertilities from brown trout in Bellbekken. 
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survival estimates (Fig. 4B). The N eV estimates were also 
rather insensitive to the shape of the age-specific proba- 
bility of maturation curve (Fig. SI), giving slightly higher 
estimate when using the age-determination-error unad- 
justed maturation probabilities (128.1 ± 51.2). 

Discussion 

In this study of an exhaustively sampled stream-resident 
brown trout population, we were able to estimate effec- 
tive number of spawners per season and generation and 
to quantify several factors that impacted the N^/N and 
N e /N ratios. This is the first study to our knowledge that 
uses parentage genetic assignments data to estimate effec- 
tive population sizes in a natural salmonid population, 
yielding separate estimates for per season and per genera- 
tion. Earlier attempts on this species have relied on 
temporal change in allele frequencies, yielding only per- 
generation estimates (Jorde and Ryman 1996; Laikre et al. 
1998, 2002; Hansen et al. 2002; Palm et al. 2003; Hegg- 
enes et al. 2009). While these approaches to estimate per- 
generation N e yielded congruent results, our approach has 
the great advantage of also providing estimates of per sea- 
son N b and of the factors determining effective sizes. As 
expected, the effect of the variance in reproductive success 
led to the largest reduction in N e . relative to the census 
size (N), whereas unequal sex ratios had only a marginal 
effect. The impact of family-size variation might be buf- 
fered with age as juvenile survival seems to be rather ran- 



dom at older juvenile stages with respect to families. 
Simulation results indicate that iteroparity, on the other 
hand, increases the variance in reproductive success, lead- 
ing to a decreased effective population size. Iteroparity 
however also decreased the variance of the effective popu- 
lation size, which might lead to a decreased risk of demise 
because of stochastic population events. 

The genetic parentage assignments and the exhaustive 
sampling should allow for robust estimates of N h . Never- 
theless, nearly half of the offspring were not assigned a 
parent, and there was some uncertainty in the ageing of 
mature fish, which might bias the generational N e esti- 
mate. We do not expect that the high proportion of 
unsampled parents would have significant impact on the 
a 2 m and of, as there is no reason to suspect that the un- 
sampled spawners produce a different family-size distri- 
bution. Also, the N e estimate seemed to be insensitive to 
ageing errors as analyses done without making age adjust- 
ments (see Materials and methods section) yielded esti- 
mates that were only slightly higher. The N e estimate also 
seems to be quite insensitive to higher family-size variance 
and survival values than the observed ones, but increase 
significantly at lower variance values (Fig. 4). This result 
is expected for the family-size variance, whereas decrease 
in survival seems to increase family-size variance and thus 
to lower generational N e value. 

We arrive at N^/Ni estimates of 0.16-0.28, whereas the 
N e /N value was around 0.5 from our simulation model. 
Because census size is often the only available demographic 
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parameter from wild populations, obtaining taxon-specific 
N h /Ni and NJN would be useful for monitoring changes 
in genetic diversity. The N b /M values are within the range 
of those reported for other salmonids, both semelparous 
(0.1-0.4, (Hendry and Stearns 2004)) and iteroparous 
(0.11-0.31, (McLean et al. 2007)). Generational N e /N is 
theoretically expected to be around 0.5 for iteroparous spe- 
cies with overlapping generations (Nunney 1993) and 
range between 0.25 and 0.75 (Nunney and Elam 1994; Wa- 
ite and Parker 1996). In a recent review, Palstra and Ruzz- 
ante (2008) reported a median N e /N value of 0.14 for a 
range of organisms, but that this ratio is actually higher for 
small populations. The formula for effective size of the 
population assumes no environmental stochasticity, stable 
age structure and constant population size. However, in 
organisms with overlapping generations, N e /N ratios below 
0.1 apparently require extreme conditions, one of which 
might be selection (Nunney 1993, 1996; Nunney and Elam 
1994). Comparisons of estimated Nb (39.8, on average) 
and N e (104) verifies that for iteroparous species N e per 
generation is not simply a product of per season N b and 
generation length (approximately 6 years), as is the case 
for semelparous species (Waples 1990). 

Family-specific survival 

Survival at later juvenile stages (older than 1+) was rather 
random with respect to families and occasionally led to 
increases in N h /Ni (Table 2). Similar pattern has also been 
found in other salmonid studies (Hedrick et al. 2000). 
More pronounced family-correlated survival might be 
expected for younger juvenile stages than what we can 
make inference for, when the juveniles belonging to a sib- 
ship group are especially vulnerable and tend to be 
spatially clumped (Serbezov et al. 2010). We observed, 
however, marked year-to-year variation in the amount by 
which N b /N is reduced between 1+ and 3+ life stages. 
Stochastic environmental or demographic events would 
explain such year-to-year variation. Year-to-year variabil- 
ity in juvenile survival has been shown to dominate over 
size-dependent survival in this system (Carlson et al. 
2008). 

Our simulation results show that the iteroparous sys- 
tem has higher variance in reproductive success that con- 
sequently lowered the N e value. This was despite the fact 
that iteroparous individuals are more productive because 
of higher cumulative fecundity and lifetime fitness (see 
also Fleming and Reynolds 2004). In polygamous systems, 
iteroparity can be thought to lead to subordinate individ- 
uals having more opportunities to be successful, reducing 
the between-individual variance in reproductive success 
and thus increasing N e (Nunney 1993). We could not 
find any significant cost of reproduction, but our results 



revealed that the effect of body size on siring offspring is 
significant (see also Serbezov et al. 2010). Larger individu- 
als leave more offspring compared to smaller individuals, 
and this difference widens as individuals reproduce 
repeatedly. Such a scenario conforms with the fact that 
family-size variance increases with decrease in survival in 
our sensitivity analysis. Other studies have also conjec- 
tured that iteroparity and multiple paternity can increase 
the variance in reproductive success and therefore reduce 
N e (Karl 2008). In the closely related Atlantic salmon 
(Salmo salar), repeat spawners can also produce dispro- 
portionate share of recruitment (Mills 1989). 

Small populations may fluctuate in numbers for a vari- 
ety of reasons, including demographic stochasticity, in 
addition to natural or human-driven environmental per- 
turbations. If population fluctuations are sufficiently large 
and not rapidly attenuated, they could lead to loss of 
genetic diversity and possibly population extinction 
(McCauley 1991). Complexity, roughly defined as the 
degree of iteroparity, is positively correlated with stability 
in terms of ability to recover from perturbations (Demet- 
rius et al. 2004). Here, the N e estimates obtained by 
assuming strict semelparity led to increased variance in 
the N e estimates, implying that such a population might 
be more exposed to stochastic demographic processes. 
There were signs of demographic/genetic stochasticity/ 
perturbations in our system. The 2005 cohort appears to 
be relatively much reduced in samplings from fall 2006, 
despite similar sampling effort and catchability (data not 
shown). The estimated total number of individuals in the 
stream varied up to 45% between the seasons investi- 
gated. This considered, and in the light of the size of the 
population, demographic stochasticity might not be 
uncommon in this system. Nevertheless, the system seems 
to be resilient to perturbations as it shows relatively high 
levels of genetic variation (data not shown). Iteroparity, 
which can be viewed as a 'bet-hedging' strategy to deal 
with environmental uncertainty, thus might be an effec- 
tive strategy in such a system which might regularly expe- 
rience years with poor recruitment and low number of 
breeders (Gaggiotti and Vetter 1999). Some support for 
this comes from empirical studies on long-lived fish spe- 
cies (e.g., Diaz et al. 2000; Lippe et al. 2006), which 
reported the maintenance of high genetic diversity despite 
relatively small population sizes. Patterns of iteroparity 
between steelhead trout (Oncorhynchus mykiss) popula- 
tions reveal that the levels of iteroparity are highest at lat- 
itudinal extremes and suggest that it might be a more 
beneficial strategy at harsher/less predictable conditions 
(Stearns 1976; Roff 2002; Einum and Fleming 2007). 

From a conservation or management standpoints, it is 
important to understand how populations might evolve 
in response to environmental changes and how these 
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responses might then feed back on population persistence 
and productivity. Among other things, this requires an 
understanding of the factors affecting effective population 
size. Here, our results emphasize the importance of 
obtaining knowledge of important life history traits and 
delineating their impact on evolutionary dynamic pro- 
cesses such as those reflected by N e . The life history and 
demographic parameters that are most essential to a popu- 
lation's conservation can thus be better identified, espe- 
cially in species with complex life histories, allowing 
researchers to design improved management and conserva- 
tion programmes. Namely, we documented large effect of 
the variance in reproductive success on reducing the N\/Ni 
ratio in a wild salmonid population. However, this 
reduction may be buffered with age of the juveniles as fam- 
ily-specific survival seems to be rather random at older 
juvenile stages. Simulation results also indicated that itero- 
parity increases the variance in reproductive success, lead- 
ing to a decreased effective population size. At the same 
time, iteroparity decreases the variance of the effective 
population size, which might lead to a decreased risk of 
demise because of stochastic population events. The evolu- 
tionary implications of effective population size therefore 
depend on detailed knowledge of its interactions with life 
history and population dynamics. These insights may be 
attained only through carefully considering age structure 
in empirical investigations of N e (Waples 2010). Finally, 
environmental change might produce conflicting selection 
pressures in different life stages in organism with complex 
life histories, which will interact with plastic (i.e., nonge- 
netic) changes in various ways (Crozier et al. 2008). In the 
case of the resident brown trout population studied here, 
we infer high levels of resilience, which results from a com- 
bination of plastic life history traits such as promiscuous 
mating system and iteroparity. Specifically, the N e seems to 
be pretty resilient with respect to age-specific maturation 
probabilities and survival rates. Nevertheless, the factors 
identified in this paper to reduce the N e should be carefully 
considered in conservation and management applications. 
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